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Abstract 

We calculate the decay rate of bottomonium to two-charm quark jets T — > cc at the tree level 
and one-loop level including color-singlet and color-octet bb annihilations. We find that the short 
distance coefficient of the color-octet piece is much larger than the color-singlet piece, and that 
the QCD correction will change the endpoint behavior of the charm quark jet. The color-singlet 
piece is strongly affected by the one- loop QCD correction. In contrast, the QCD correction to 
the color-octet piece is weak. Once the experiment can measure the branching ratio and energy 
distribution of the two-charm quark jets in the T decay, the result can be used to test the color 
octet mechanism or give a strong constraint on the color-octet matrix elements. 

PACS numbers: 12.38.-t, 12.39.St, 13.20. Gd,14.40.Gx 
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I. INTRODUCTION 



It is commonly believed that the heavy quark pair production and annihilation decay can 
be described by perturbative Quantum Chromodynamics (pQCD) since the heavy quarko- 
nium mass provides a scale that is much larger than A.qcd- Due to its nonrelativistic nature, 
the heavy quarkonium annihilation decay is expected to be described in an effective theory, 
non-relativistic Quantum Chromodynamics (NRQCD) [lj . In the NRQCD factorization for- 
malism, the decay of heavy quarkonium is described by a series of annihilations of the heavy 
quark pair states and corresponding long-distance matrix elements, which are scaled by the 
relative velocity v of quark and antiquark in the quarkonium rest frame. The heavy quark 
pair states can have not only the same quantum numbers as those of the quarkonium, but 
also other different quantum numbers in color and angular momentum. In particular, the 
heavy quark pair can be in a color-octet state. 

The color-octet scenario seems to acquire some significant successes in describing heavy 
quarkonium decay and production. But recently, several next-to-leading order (NLO) QCD 
corrections for the inclusive and exclusive heavy quarkonium production in the color-singlet 
piece are found to be large and significantly relieve the conflicts between the color-singlet 
model predictions and experiments. It may imply, though inconclusively, that the color- 
octet contributions in the production processes are not as big as previously expected, and 
the color-octet mechanism should be studied more carefully. 

The current experimental results on inelastic Jlip photoproduction at HERA are ade- 
quately described by the NLO color singlet piece [2J . The DELPHI data favor the NRQCD 
formalism for J/ip production in 77 — > J/ipX, rather than the color-singlet modelfl, H[. 
The large discrepancies in J/ip production via double cc in e + e~ annihilation at B factories 



are 



between LO theoretical predictions [5], @, 0, S, @, 13] and experimental results 11, 12 
probably resolved by including the higher order corrections: NLO QCD corrections and 
relativistic corrections [l3|, 14, 15, lfj 13, 18, 19]. The NLO QCD corrections in J/ip and T 



production at the Tevatron and LHC are calculated including the color singlet piece [20j, [21 
and the color octet piece [22{|. The QCD corrections to polarizations of J/ip and T at the 



24J . The experimental data of polarizations 



Tevatron and LHC are also calculated |22|, [ 
at the Tevatron seem to favor the NLO QCD corrections of the color singlet piece rather 
than the color octet piece. Recent developments and related topics in quarkonium physics 



can be found in Refs. 25, 26, 27 



In order to further test the color octet mechanism, in this paper we calculate the rate 
of bottomonium decay into a charm quark pair T — > cc. There have been some works 
on bottomonium decays and the color octet mechanism. Fritzsch and Streng calculated 
the decay rate of T into charm at leading order in a s , T — > ggg* — > ggcc\28]. Bigi and 
Nussinov have taken into account the contribution of T —>■ gg*g* — ► gcc [29j. Barbieri, 
Caffo, and Remiddi have calculated the decay rates of the P-wave bottomonium states 
into charm at leading order in a s [30] . Maltoni and Petrelli calculated the effects of color- 
octet contributions on the radiative T decay [3l[. Recently, Bodwin, Braaten and Kang 



calculated the inclusive decay rate of Xb into charmed hadron in the NRQCD framework [32 
Gao, Zhang and Chao calculated the bottomonium radiative decays to charmonium and 
light mesons|33, |34J, as well as T radiative decay to light quark jet to test the color octet 



mechanism I35I1 
and a? I36l, \3j 



The S-wave quarkonium decay to light hadrons was calculated up to order v 4 
The exclusive double charmonium production from T decay was calculated 



by Jia [38j. Kang, Kim, Lee and Yu have calculated the inclusive charm production in T(nS) 
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decay |39J]. And the invariant-mass distribution of cc in T(l S) — > cc + X was also calculated 



by Chung, Kim and Lee|40j. And r]b inclusive charm decay was calculated by Hao, Qiao 
and Sun 41|. As to experiment, ARGUS Collaboration searched for charm production in 
direct decays of the T(1S), and found B dir [T(lS) -> D*(2010) ± + X] < 0.0190. Very 
recentl y C LEO has searched for the D° production in direct decays of the Xbj{ n S) (n=l,2) 
states [43J. The present investigation for the T decay to cc pair will hopefully add a new 
contribution to the test of color-octet mechanism in heavy quarkonium decays. 

This paper is organized as follows. In Sec. [Til we present the theoretical framework for 
the decay of T — > cc. In Sec. IHh we estimate the color-singlet contributions. In Sec. [TV] we 
include the color-octet contributions. In Sec.|V] we discuss the NRQCD matrix elements e.g. 
(T|(9( 3 5'i > 8) |T) and give a numerical estimation of the color-octet contributions. Summary 
and discussion are presented in Sec. [V]] The detailed and lengthy intermediate steps and 
formulas in the calculation will be given in the appendices. 



II. THEORETICAL FRAMEWORK 



In the framework of NRQCD, the width of T decay to cc can be written as 
T[T -> cc] =^f[66(n) -> cc] (T|0(n)|T), 



(1) 



where n denote quantum numbers including the spin angular momentum S, orbit angular 
momentum L, total angular momentum J, and the color index 1 or 8. The short distance co- 
efficients T[bb(n) — > cc] can be calculated in pQCD, and the long distance factors (T|C(n)|T) 
scale as definite powers of the relative velocity v of quark and antiquark in the quarkonium 
rest frame For T, the leading order matrix element is (T|(9( 3 S'i j i)|T), and there are 
three matrix elements that contribute up to corrections of relative order v A \ (T|C( 1 S , o i s)|T), 
(T|(9( 3 5'i i 8)|T), and (T|(9( 3 Pj 8 )|T). The other matrix elements are of higher order in v. 




bb( 3 S 1A 




FIG. 1: Feynman diagrams for the color-singlet decay T(66( 3 S'i ) i)) 
and three virtual gluons (right). 



cc via a virtual photon (left) 



Feynman diagrams for the color-singlet decay T(bb( 3 Si^)) — > cc via a virtual photon (left) 
and three virtual gluons (right) are shown in Fig. [TJ At leading order in a s , the decay of the 
color-singlet piece T(66( 3 S'i i i)) can proceed through a virtual photon or three virtual gluons. 
The decay width is of order 0((a/n) 2 ) for the virtual photon, and 0((a s /n) 6 ) for the three- 
gluons. So the single photon process is expected to be dominant, and the contribution of 
three-gluon process will be roughly estimated in Sec. IIHI If a soft gluon is allowed to appear 
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bb( 3 S 1>8 , 

5*0,8; Pj, 




FIG. 2: Feynman diagrams for the color-octet bb — ► cc 



in the final state, the order of a s in the process can be decreased. But the processes of 



bbrs 



2g* + g — ► cc + g and bb( £1,1) — >■ g* + 2g — >■ cc + 2g are infrared (IR) finite 



28l . [29|, so the phase space of the soft gluon will bring a suppression factor: 



d ka 



m l k g 



(2) 



k° g <m b 5s 



where the factor of ml is used to balance the dimension, and 5 S is the soft cut. The gluon 
is regarded as a soft gluon when the energy of the gluon is lower than rrib5 s . When 5 S is set 
to, say, 0.2, the corresponding energy cut is about 1 GeV, 5 2 s is numerically close to a s /n, 
so these soft gluon processes are relatively suppressed and should be ignored here. 

The decay of the color-octet piece bb — > cc includes contributions from color-octet bb 
components 3 Si )8 , as well as 1 S 0j $ and 3 Pj > 8 in the T Fock state expansion. Feynman diagrams 
for the color-octet bb — > cc are shown in Fig.[2j The leading order decay width of bb( 3 Si t s) 
cc is of order 0(a^/7i 2 ), while processes 66( 1 i So,8, 3 Pj,s) cc can only proceed via a loop, and 
the corresponding decay widths are of order 0(a^/ir 4 ). Moreover, since (T|(9( 3 S'i ) 8)|T) ~ 



<T|0(% i8 )|T> 



(T|o( 3 P Ji8 )jr) 



(T|C( 3 S'i i i)|T) according to the velocity scaling rule, 



the contributions of T(66( 1 S'o i 8, 3 Pj,s)) cc can be neglected. 

The color-singlet and color-octet contributions will be discussed respectively in the next 
two sections. 



III. COLOR-SINGLET PIECE bb( 3 S 1A ) -> cc 

The amplitude of color-singlet piece bb( 3 Si t i) — > cc can be written as [33|, 135 



A[bb( 3 S lil (2 Pb )^c( Pc ) + c(p- c 



(T\0(?S 1A )\T) 



2N r 



E EE 

LrzSrz s l s 2 jk 



x(l I 3k;3j)(J r J Tz I L r L Tz ;S T S Tz )(S r S rz \ s 1 ;s 2 ) 
xA(bj(p b ) + b k {p b ) -> ciip c ) + c~i(p 5 



(3) 



where (1 | 3k; 3j) = 5jk/\fN~ci (SrSxz I si; s 2 ), and (JrJrz \ L^L^ Z ; SrS^z) are respectively 
the color-SU(3), spin-SU(2), and angular momentum Clebsch-Gordan coefficients for bb pairs 
projecting on appropriate bound states T. A(bj(pb)+b~k(Pb) ~* c j(Pc)+Ci(Pc)) is the am plit ude 
of the process bj(pb) + bk(Pb) q(Pc) + 5i(Pc)- I n the calculation, we use FeynArts |44j, |45 
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to generate Feynman diagrams and amplitudes, FeynCalc 46j for the tensor reduction, and 
LoopTools [47] for the numerical evaluation of the IR-safe integrals. 

The spin projection operators PssAVil) which describe quarkonium production are ex- 



pressed in terms of quark and anti-quark spinors as[48l. 149 



Sl,S 2 



Pss z (p, q) =^y{^+q,s 2 )v^-q,s l )(s 1 ;s 2 \S S z ) 



For the 3 5 , i state, it is 



= (2 A + 2m) jt{S z ). 



(4) 



(5) 



And the spin projection operators which describe the annihilation of quarkonium are the 
complex conjugate of the corresponding operators for production. 

The leading order (LO) color-singlet decay bb( 3 Si t i) — > 7* — ^ cc is shown in Fig. [TJ and 
the corresponding LO decay width is 



r i0 [T( 3 S' 1 , 1 ) 



cc 



47ra 2 v / T^(2 + r 2 ) 



TlOfSVOlT) 



(6) 



where r 2 = m 2 /m 2 . This result is consistent with Ref. 39|, |40j. Comparing it with the 
leptonic width 



r i0 [T 



e e 



2na 2 
27m 2 



;t|0( 3 s , 1j1 )|t>, 



(?) 



we can get 



r L o[T( 3 5 1 , 1 ) - cc] = - r LO [T - e+e"] X (l + O (r 2 )) 



(8) 



where the factor 4/3 comes from the charm quark charge and color factor, and r 2 
If we set m& = 4.7 GeV and m c = 1.5 GeV, the LO decay branching ratio is 



10" 



B LO [TCS 1A 



cc\ 



1.33S i0 [T -> e+e"] = 3.: 



(9) 



where B[T -> e+e"] = (2.38 ± 0.11)% is used according to the PDG 2006 Version|50j. 

We next consider the QCD radiative corrections. The Feynman diagrams of one-loop 
virtual corrections and counter terms are shown in Fig. [3j The renormalization of heavy 
quark wave function should appear. The on-mass-shcll (OS) scheme is chosen for Z 2 b and 

z 2 M 



5Z? S 



J 2b 



5Z, 



OS 
2c 
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euv em ra% 
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+ 0(a 



2 ) 



(10) 



where \x is the renormalization scale, je is the Euler's constant. In this scheme, we need not 
calculate the correction of external quark legs. We employ the two-loop formula for a s (/i), 



fa InL 



4vr (3 L /3 3 L 2 ' 



(11) 
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Real Nl Real N2 



FIG. 3: Feynman diagrams for one-loop QCD corrections with counter terms for 66( 3 S'i i i) — » 7* — ► 
cc 

where L — In (fi 2 /A 2 QCD ), and ft = (34/3)C A 2 - AC F T F n f - (20/3)C A T F n f is the two-loop 
coefficient of the QCD beta function. 

The correction to T(bb( 3 Si^)) — > 7* gives a factor of — l|p>. a t 0(a: s ). The other part is 
the correction to 7* — > cc. If we set m c = 0, then the combined total correction becomes the 
correction to R, the ratio of cross section of e + e~ — > Zig/it hadrons to that of e + e~ — > ^ . 
It give a factor of — at 0(a a ). 

Compared with it, the leptonic width of T at 0(a s ) becomes the known result: 

rWT - e + c-] = |^ (l " ^) <T|O( 3 Si,0|T) . (12) 

If the parameters are chosen as mt, = 4.7 GeV, m c = 1.5 GeV, and a s = 0.220, then the 
branching ratio is 

B NLO [T( 3 S ltl ) ^cc + X] = l.QB NLO [T -> e+e^] w 3.8%. (13) 

When the emitted gluon energy is large, it would form a jet. So a cut to the gluon energy 
should be introduced to distinguish between cc and ccg final states. If E g < nib x 5 S , the 
gluon is considered as soft and the final state is cc. Otherwise, when E g > x 8 S , the 
final state is ccg. If we set 5 S = 0.15, then the branching ratio is 

iWT^Si,!) - cc) w 1.4%. (14) 
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If we set 5 S = 0.10 and 0.20, then the branching ratio is about 5.4 x 10~ 3 and 2.0% respec- 
tively. 

For the color-singlet piece, the contribution of T( 3 S'i ) i) — > 3g* — > cc, which is shown on 
the right hand side in Fig. [TJ has not been calculated so far, and we may have a rough 



estimate for it. We can use T( s Si 



3g to give an order of magnitude estimate for the 



contribution of Y( 3 Si ; i) — > 3g* — > cc. We have the following order of magnitude estimates: 



B[T( 3 S 1A ] 



B[T( 3 £i,i) 
) ^ 999 -» 



S[T( 3 5 ltl 



/+/-] oc 

-> 3flr] oc 

#5>cc] oc 

— > ccl oc 



7T 



Comparing the leptonic width with three-gluon width, we get 



rprgj^o -> /+/-] (f) 2 
r[T(3^ M ) - 3^] (^) : 



0.016, 



which is about one half of the experimental value of r[T( 3 5i i)->3g] ~ 0-03 ^Oj, anc ^ ^ e 
closeness of this estimate to the data may suggest that the naive estimate could make sense. 
Comparing T( 3 S'i i i) — > 3g* — > cc with T( 3 S'i ] i) — > 3g and T( 3 5'i i i) — > 7* 



(15) 



(16) 



r[T( 3 s 



1,1/ 



3g* 



cc\ 



r[r(*s hl ) - 3^] 



^3x 10- 4 , 

7T / 



cc, we can get 
(17) 



r[T( 3 5 1>1 ) -> 3£* - cc] c^ 
T[T( 3 5 1 , 1 ) -> 7* - cc] 



a 2 7r 4 



0.02. 



(18) 



From the estimates given in Eq.( fT7|) and Eq.t flBl . we see that the contribution of T( 3 S'i ) i) — > 
3g* — > cc is very small and much smaller than that of T( 3 Sx,i) —>■ 7* — ► cc. Even if the 
contribution of T( 3 S'i ) i) — > 3g* — ► cc is underestimated by an order of magnitude in Eq. (fT7|) 
and Eq. (fT8|) . we could still expect that for the decay T( 3 5i t i) — > cc the QED contribution is 
dominant. Another useful example is the decay rate rjTf^S'-i.i ) — > ggcc], which is of higher 
order in a s than r[T( 3 S'i ) i) — > ggg], and is given in Ref . [391] . From their estimate we can get 

r^( 3 5i'o^lggg^ = 0-029, which is also of the same order of magnitude as, but even smaller 



than, our naive estimate. 



r[T( 3 s ; 



1,1/ 



ggcc] 



1,1, 



ggg] 



V 7T 



0.071. 



(19) 



This might imply that the naive estimates given in Eq.( jT5j) as well as in Eq. tflBl might 
be tenable in estimating the rates of higher order processes by order of magnitude. So, 
based on the rough estimate given in Eq. (1181) for the contributions of the color-singlet piece 
to the T( 3 5'i i i) — ► cc process, we assume that, as an approximation, the contribution of 
T( 3 5'i,i) — > 3g* — > cc can be neglected, and only T( 3 S'i ) i) — > 7* — > cc will be taken into 
consideration. 
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IV. COLOR OCTET PIECE bb( 3 S 1:8 ) cc 

The amplitude of color-octet piece &&( 3 Si i8 ) — > cc can be written as |33|, 135 



A(bb( 3 S^(2p b ) ^c( Pc ) + c(p. 



(T|0( 3 5 li8 )|T> £ EE 

LtzStz s 1 s 2 jk 

x(8 a | 3k;3j){JJ z \ LL Z ; SS Z )(SS Z | s x ; s 2 ) 
x^4( 6j(p 6 ) + 6 fe (p 6 ) -> Q(p c ) + Ci(p £ 



(20) 



where (8 a | 3k; 3j) = y/2T? k , and other expressions are similar to the color-singlet piece. 

The Born diagram of &6( 3 Si i8 ) — > cc is shown in Fig. [2j It is also calculated in Ref. 32 
The leading order width is 



^' ^ 6^ +r> (T| ° (3g '' 8)|T> - 



(21) 



We further calculate the next-to-leading order (NLO) corrections. The Feynman diagrams 
for NLO virtual corrections with counter terms in the color-octet piece 66( 3 5i. 8 ) — > cc are 
shown in Fig. HI The Feynman diagrams for NLO real corrections in the color-octet piece 
66 ( 3 5*1,8) — ► cc are shown in Fig. [5j The renormalization of heavy quark wave function, gluon 
wave function, and coupling constant should appear here. Z 2 & and Z 2c are given in Eq. ffTOl) . 
For Z 3 and Z g , we choose the modified minimal-subtraction ( MS ) scheme 17 



5Zf s 



a 



^(/?o-2CU) 



5Z 



MS _ 

2 4tt 



7s + ln(47r) 



7s + 1h(4tt) 



2 ) 

s ) • 



(22) 



The parameters are chosen as m b = 4.7 GeV, m c = 1.5 GeV, rif = 4, Aq^ d = 338 MeV, 
/i = mfe, and then a s = 0.220. So we can get the leading order result: 



5lo[T( 3 5 1i8 ) 



cc 



42 x 



(T|0( 3 5 l!8 )|T) 



The total NLO result is 



-6jVLo[Y( 3 5i,. 



cc + X] = 53 x 



GeV 3 

(T|Q( 3 g 1|8 )|T) 
GeV 3 



If we set the soft cut 5 S = 0.15, then the NLO result is 

(T|G( 3 5 li8 )|T) 



Bnlo [^( 3 5i, 8 ) 



cc 



41 x 



GeV 3 



(23) 
(24) 

(25) 
and 



/ y. | ^) ^3 i^j.^ g ^ | 

If we set the soft cut 5 S = 0.10 and 0.20, then the branching ratio is 37 x j — — - 

.. (t|0( 3 5i, 8 )|t) 

^4 x — GeV 3 — respectively. 

From the above expressions, we see that the short- distance coefficient for this color-octet 
process is large, and this color-octet process may make a significant contribution to the T 
decay to two-charm quark jet. The numerical estimate will be given in the next section. 
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CT Nl 



CT N2 




Box Nl Box N2 

FIG. 4: Feynman diagrams for next-to- leading order virtual corrections with counter terms in the 
color-octet piece bb( 3 S\ t s) — > cc 



The color-octet pieces 66( 3 Pj j8 ) and bb^Sog) also contribute to the charm quark jet 
production through ccg, where the gluon is soft. The bb( 3 Pj^) —>■ ccg is IR divergent, and 
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bb( 3 S_ 




Real N5 

FIG. 5: Feynman diagrams for next-to-leading order real corrections in the color-octet piece 
bb( 3 S 1;8 ) -» cc 



it should be absorbed into the matrix element (T|0( 3 S , i i g) |T) 



[t\o( 3 s 1>8 )\t) 1 = (r\o H ( 3 s 1 ,s)\r) 



Ca\ TTOs 

2 ) 2v 



+ 



4a s / 47r/i 



,2\ e 



A 2 



x 



exp(-e lE ) (-L-±-)Y2 Bf <T |C ( 3 P,, 8 ) | T> , 
V^uv e m J j=q 



(26) 



where the Coulomb term of (T |(9 (^i^)! T) 1 is canceled by the virtual correction, the IR 
divergent terms is canceled by bb( 3 Pj 8 ) — ► ccg, and the UV divergent term gives the running 
of matrix element. If we choose the matrix element renormalization scale as m&, then we 
find the branching ratio of bb( 3 Pj 8 ) and bb^So^) decays into ccg at order of a 3 to be 



B[rCs , 8 

P[T( 3 P 



cc + X] 



2.8 x 



J,8) 



cc + X] = 0.61 x 



GeV 3 
(T|0( 3 P o , 8 )|T) 

GeV 5 



(27) 



Since (T 1 0(^0,8) | T), (T|0( 3 Si, 8 )|T), and (X |0( 3 Pj, 8 ) |T) jm\ are of the same order , we can 
ignore the contribution of fefe( 3 P/ j8 ) and bb^So^), as compared with the bb( 3 Si^) contribution 
given in Eq. fT24|) . 
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FIG. 6: Decay widths of the color-singlet piece and color-octet piece bb( 3 Si) — > cc rescaled by 
the corresponding value at m c = 1.5 GeV as functions of the charm quark mass m c . Here A = 
0.338GeV, m& = 4.7 GeV, \i = m^, and the soft cut 5 S = 0.15. LO means leading order, and NLO 
means next-to-leading order. 3 S'i i i means the ratio T[bb( 3 Si t i — > cc](m c ) /T[6fe( 3 Si t i — > cc](m c = 
1.5 GeV) in the color-singlet piece, and 3 5'i i 8 means the corresponding ratio in the color-octet 
piece. 



The dependence of the leading order and next-to-leading order decay widths in the color- 
singlet and color-octet pieces bb —>■ cc on the charm quark is shown in Fig. [61 The dependence 
of LO result on the charm quark mass is weak and the same for the color-singlet and 
color-octet pieces. The reason can be found in Eq. (jSJ) and Eq. ( 12"T|) . If we choose m c = 
1.5 ± 0.2 GeV, the ratio is about 1 ± 0.003 at LO in a s , lloJo a ^ NLO for the color-singlet 
piece, and llo!o57 a ^ NLO for the color-octet piece. 

The dependence of the leading order and next-to-leading order decay widths in the color- 
singlet and color-octet pieces bb —>■ cc on the renormalization scale fi is shown in Fig. [71 The 
LO color-singlet result is independent of the renormalization scale. As it is shown in the 
curve of NLO 3 5'i i 8, we choose \x = vn^ for the principle of minimum sensitivity(PMS) [5l| . 

The dependence of the next-to-leading order decay widths in the color-singlet and color- 
octet pieces bb — > cc on the soft cut 5 S is shown in Fig. [HI The LO result is independent 
of the soft cut 5 S . The NLO color-singlet result is rather sensitive to 8 S , whereas the NLO 
color-octet result is insensitive to 5*. 



11 



NLO ^i,! 

LO 3 5i,i 
NLO 3 S 1;8 

LO 3 5 1)8 

fi(GeV) 

1111111111111111111111111111111111111111111111 
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FIG. 7: Decay widths of the color-singlet piece and color-octet piece bb( 3 Si) — ► cc rescaled by the 
corresponding value at /i = m& as functions of the renormalization scale ^. Here A = 0.338GeV, 
nib = 4.7 GeV, m c = 1.5 GeV, and the soft cut 5 S = 0.15. LO means leading order, and NLO means 
next-to-leading order. 3 61,1 means the ratio of r[66( 3 5i 5 i — > cc](//)/r[66( 3 S'i i i — > cc](/i = m&) in 
the color-singlet piece, and 3 5'i i 8 means the corresponding ratio in the color-octet piece. 




V. COLOR OCTET MATRIX ELEMENTS 

The color-singlet matrix element (T|(9( 3 5'i i i)|T) can be extracted from the T leptonic 
decay width. Using Eq. f|T2l) . we get 

(T\0( 3 S ltl )\T) = 3.8 GeV 3 . (28) 

On the other hand, large uncertainty is related to the color octet matrix element 
(T|C( 3 S'i i8 ) |T). According to the velocity scaling rule and taking v 2 = 0.08, we might 
naively have 

<T|G( 3 5 li8 )|T) w ^ <T|0( 3 5 1 , 1 )|T) = 4.1 x HT 3 GeV 3 . (29) 

Using Eq. (j2"4"|) . we would get 

B NLO [T( 3 S h8 ) -> cc + X] = 21%. (30) 

For the light quark q = u,d,s, we have 

Stvlo[T( 3 5 1i8 ) ^ qq + X] = B NLO [T( 3 S ltli ) - cc + X] x (l + (r 2 )) . (31) 

So T decay through bb( 3 Si.$) would have a very large branching ratio, say about 80%. 
Apparently, the color -octet matrix element estimated in this naive way from the velocity 
scaling rule is greatly overestimated, even by an order of magnitude. 
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CM 



NLO 3 S h8 
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0.4 
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FIG. 8: Decay widths of the color-singlet piece and color-octet piece bb( 3 Si) — ► cc rescaled by 
the corresponding value at 5 S = 0.15 as functions of the soft cut S s . Here A = 0.338GeV, nib = 
4.7 GeV, \x = nib, an d ni c = 1.5 GeV. NLO means next-to-leading order. 3 <Si i means the ratio of 
r[66( 3 S'i > i — ► cc](5 s )/r[66( 3 S'i j i — > cc]((5 s = 0.15), and 3 £1,8 means the corresponding ratio in the 
color-octet piece. 



Another approach to determine the matrix element is the lattice QCD calculations. The 



lattice calculation in Ref. 52j gives 

(T|0( 3 £i )8 )|T) « 8.1 x 10" 5 <T|C( 3 Si,i)|T) = 3.1 x 10~ 4 GeV 3 . (32) 

If we set the soft cut 5 S = 0.15, then the next-to-leading order result is 

B NLO [T{ 3 S 1>8 )^cc] = 1.3%. (33) 

If we set the soft cut 5 S = 0.10 and 5 S = 0.20, the branching ratio is 1.1% and 1.4% 
respectively. 

From the above numerical results and Eq. (1241) we see that since the short distance coeffi- 
cient for the color-octet contribution to the T — > cc decay is large, this process is sensitive to 
the value of the color-octet matrix element, and may therefore serve as a useful test ground 
of the color-octet mechanism. 

Moreover, the next-to-leading order QCD correction in the color-singlet piece is much 
stronger than that in the color-octet piece, and the color-singlet contribution shows a strong 
sensitivity to the soft cut 8 S , whereas the color-octet result does not. These differences 
between the color-singlet and color-octet contributions will also be significant in clarifying 
the issue about the color-octet mechanism. Once the experiment can measure the branching 
ratio and energy distribution of the charm quark jet in the T decay, the result can be used to 
test the color-octet mechanism or give a strong constraint on the color-octet matrix elements. 
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VI. SUMMARY AND DISCUSSION 



We calculate the decay rate of bottomonium to two-charm quark jets T — > cc at the tree 
level and one-loop level including color-singlet and color-octet bb annihilations. We find that 
the short distance coefficient of the color-octet piece is much larger than the color-singlet 
piece, and that the QCD correction will change the endpoint behavior of the charm quark 
jet. The color-singlet piece is strongly affected by the one-loop QCD correction. In contrast, 
the QCD correction to the color-octet piece is weak. Once the experiment can measure 
the branching ratio and energy distribution of the two-charm quark jets in the T decay, 
the result can be used to test the color octet mechanism or give a strong constraint on the 
color-octet matrix elements. 



After our work was completed [53|, a paper appeared [39J, in which Kang, Kim, Lee, 
and Yu calculated the inclusive charm production in T(nS) decay. They focused on the 
inclusive charm production of the color-singlet piece at leading order in the strong coupling 
constant a s . We focused on the cc final state and the color-octet mechanism. The cc final 
state is essentially the two charm-jet process. And we have calculated the next-to-leading 
order QCD corrections in both color singlet and color octet pieces. Our leading order result 



of T — * 7* —> cc is consistent with their result [39 
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APPENDICES 



APPENDIX A: THE SCALAR FUNCTIONS 

The scalar functions that appear in the virtual corrections are listed in this Appendix. 
There are UV, IR and Coulomb singularities in the scalar functions. The UV and IR singu- 
larities are regularized with D = 4 — 2e space-time dimension. The exchange of longitudinal 
gluon between massive quarks in Vertex Nl in FigJ3]and FigJUleads to a Coulomb singularity 
~ 7r 2 /f, where v = \/—{pb — Vl) 2 l m b is the relative velocity between b and b in the meson 
rest frame (v — \pb — pt\/rrib). The Coulomb singularities should be canceled by that in the 
matrix elements (see, e.g., flq. fl~7l]). 

Since the imaginary part of the integrals will disappear in the final result, only the 
real parts are given. The external particles are taken to be on-mass-shell, pi 



Pc 



m c Pb-Pc= Pb-Pc 



m 2 ,, and p c ■ p s 
The scalar one-point function is defined as 



2m 2 



Pi 



77 i 



b- 



A (m 2 



A-D 



d D q 



(2tt) d q 2 



m" 



iC e {m) m 2 



1 

- + 1 

e 



(Al) 
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where 



C e (m) 



L_ e -e(7s-ln47r) f ^ 



16tt 2 



2 \ e 



(A2) 



and D = A — 2e. 

The scalar two-point function is defined as 

d D q 



B (p,m ,m 1 ) = fi 4 D J 



1 



(2n) D [q 2 - ml] [(q + p) 2 - m\] (A3) 
Four different types of two-point functions appear in the calculation of the virtual corrections: 



B (p b ,0,m b ) = B (2p b ,m b ,m b ) = iC e (m b ) 

B o (p c ,0,m c ) = B (pc, 0,m c ) = iC e (m 

B (2p b ,m c ,m c ) = iC e (m c 

Bo(p c -Pb,m c ,m b ) = iC e (m c 



i + 2 

e 



1 

- + 2 

e 



e V 1 + P 



B (2p b ,0,0) = iC e {m c 



1 n 1, (1-13 
- + 2 + - In - 

e (3 \l + (3 

1 , /4m? . 
In ( — ^ ) + 2 



Here and below we will use the shorthand notation (3 = y/l — r 2 = a/1 — m 2 / 
The scalar three-point function is defined as 

Co(pi,P2,m ,m 1 ,m 2 ) = 



mi 



A-D 



d D q 



(27r) D [q 2 - ml] [(q + pi) 2 - m\] [(q + p 2 ) 2 - m 2 ] 
The following types of three-point functions appear in the virtual corrections: 



(A4) 
(A5) 
(A6) 
(A7) 
(A8) 



(A9) 



Co(pc, -p 5 ,0,m c ,m c ) 



C (p c ,p b ,0,m c ,m b ) 



Co(Pc -Pe,m c , 0,0) 
C (Pb, -Pb, m b , 0, 0) 



iC t (m c ] 
Ami (3 



lnx/3 — 21nx i gln(l — xp) — 2Li 2 (x / g) 



+ -ln 2 x /3 -4C(2) 



iC e (y/m c m b ) 
Am c m b x/(x 2 ~ 1) 
4x 



-2 In 



iln 

e 

In 



x+iy 2 Vx+i 



(A10) 



In 2 r 



-Li 2 ^1 + 
1 

U(47r) 2 m 2 /3 
. In 2 
'(47r) 2 m 2 ' 



i-x 

(x + i) 2 ^ 11 Vx + i 

r(x- 1) 



Li 



- Li 2 1 + 



(x-i) 2 
(x + i) 2 

X-l \ vr 2 
r(x + l)J + 6 



2Li 2 (-x /3 ) + -ln 2 ^ + C(2) 



(All) 
(A12) 
(A13) 
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where x@ — (1 — + (3) , £(2) = vr 2 /6, r = m c /m b , and x = a/(1 — r)/(l + r). There is 

another scalar three-point function that is IR and Coulomb divergent, 



C (Pb, -p- b ,0,m b ,m b ) 



.C e (m b ) 
2m 2 



1 7T 2 

- + — - 2 + 0(e) 
e f 



(A14) 



where v = \/—(pb — Vl) 2 l m b- In the meson rest frame, we have v = \pl — pi\/m b . 
The scalar four-point function is defined by 

D {pi,p 2 ,P3,m , m 1 ,m 2 ,m 3 ) = 

d D q 1 



A-D 



(27r) D [q 2 - ml] [(q + pi) 2 - raf] [(g + p 2 ) 2 - m 2 ] [(q + P3) 2 - ml] 
There are three different types of four-point functions: 



(A15) 



D (p b , Pb - Pc, ~Pb, m b , 0, m c , 0) 
iC e (^m c m b ) f (x 2 - 1) 



-2 In 



4x 



-Li 2 1 + 



(x + i) 2 
x- 



rx 



In 



iln 

e 

- X 



r(x + 1) 
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+ 



7T 



6 



i-x 
x + i 

- Li 2 

1 





-In 2 



(x 



(x + i) 2 



i-x 
x + i 



- Li 2 1 



In 2 r 



2 

Kx - 1) 
x + i 



2Li 2 (-^) + -ln 2 ^ + C(2) 



(A16) 



The IR and Coulomb singularities can be regularized by the gluon mass m g . The rela- 
tion between the gluon mass m g regularization and the dimensional regularization for IR 
singularity is 



1 , 4vr/i 2 
7e + m ■ 



(A17) 



e m 

And the relation between different regularization schemes for the Coulomb singularity is 

2irm - 1 



X 



7T 
V 



(A18) 



Eq. flA17D and Eq. (IA18j) are consistent with Ref.[2j. 



APPENDIX B: REAL CORRECTIONS AND THE THREE-BODY PHASE 
SPACE 

For the real corrections, the process T(2p b ) — > c(p c ) + c(p 5 ) + g(k) is a three-body decay 
process. Similar to the method in Ref. [47J], we can write down the Lorentz-invariant phase 
space 

dPS 3 (2p b ;k, Pc ,p,) = {2 ^ k0 (2 ^ g ( ^g pg {2,f5\2 Pb -k- Pc - Pe ). (Bl) 
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Introduce the identities 

d 3 Pi 
2rf = 



d 4 p i 5(p 1 2 - mf 



|pt| 2 d\pi\ dQi \pi\ dp? df2; 



2p° 



(B2) 



where is the mass of particle i, and dfli is the direction of particle i in the 3 dimension 
space. Then we can rewrite dPS 3 (2p b ; k, p c , p g ): 



dPS, 



\k\\p c 
4(2vr) f 

I k \\p c 
4(2tt) ; 



dk° dQ g dp° c dQ c d% 5(pl - ml) 5 4 (2 Pb -k-p c -p d ) 

dk° dQ g dp c dtt c 5[{2p b -k- p c f - m 2 ] (B3) 



We define the momenta in the rest frame of the T, 

2p b = (2m b , 0,0,0) p c = (p°, \]%\ sin 0,0, cos6») 

fc = 0, 0, \~k\) Pe = (p%-\p^\wie,0,-\k\-\p^\coaO), (B4) 



where is the angular between g and c , and \pl\ — \J (p°) 2 — m 2 . Then gives a factor 
47r, and (ifi c = dcosOdcf) and <i0 gives a factor 27r. So we have 



dPS 3 = 



\_k_\\pl 
2(2ixf 



dk° dp° c dcosO S[(2p b — k — p c 



f-m 2 }. 



(B5) 



Then we use the 5 function to remove 9 in the integral with 

(2p b - k - p c ) 2 - m 2 = A;°-p°) 2 - | 2 + \p^\ 2 + 2\~k ||^|cos0) - m\ 



and 



and get 



EE /(cOS0) 

df (cost 



dcos9 



2| A; ||pc|, 



(B6) 



(B7) 



cos# 



(yfs -k° - pi) 2 - | k 



\Pc 



\ 2 — m 2 



2| A; ||p c 



and 



dPS, 



4(2tt) ; 



dk° dpi . 



(B8) 



(B9) 



To determine the limits of integration, we employ the restriction of | cos#| < 1 and p° > mf, 
then we get 



(k°) mm 
(k°) maX = m b 



a- 



(m c + m g ) 2 — m 2 
4m b 



(B10) 
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and 



WD 



max, mm 



2r 



o[j + m + mJ) ± | k U / (r — m+)(r — m_) 



r = a — \ k 



m± = m c ± nil 



(Bll) 



Here we keep the gluon mass m g for massive gluon regularization. There is a soft divergence 
in the real corrections, so we should introduce a soft cut E s for the gluon. Then the phase 
space is divided into two regions: 



dPS, 



dPSr ft | k o <Es + dPSf ard | k0>Es 



(B12) 



The hard region can be integrated in four dimension or with massless gluon. And the phase 
space in the soft region is 



(iPSf oft L p = dPS 2 / dn^- 1 I l kl °\ dk°. 

■s \k°<E s z I g / 2(2n) D ~ 1 



The decay amplitude of the color singlet process can be written as 

A Real (bl(3 Slil (2 Pb ) - C { Pc ) + c(p g ) + 9 {k) 



k°<E s 



g s ^eUk)A Born bb(*S hl (2p b ) -> c(p c ) + c(p c -) ® T a 



P^ 



Ps 



Pc ■ k Pc- k 



and 



A 



Real I 3 



<->l,l 



k°<E s 



where 



1 ■ 



Pi ■ Pj 



Pi ■ kpj ■ k 

The decay amplitude of the color octet process can be written as 

A Real (MfS^QpJ - c{ Pc ) + C{p 5 ) + g{k) 

and 



k°<E s 

P ^ -T a ® A Born - A Born ® T a P ' 



p c - A; 



Pe • k 



A 



Born 



2 p b -k 



(B13) 



(B14) 



A Born | 2 g 2^ ^ _ ^ + ? (B15) 



(B16) 



(B17) 



A 



Real ( 3 



kO<E s 



A 



Born I 3 



'Si, 



2 2e 

9s V 



4 1 4 

g-^cc "I - T^cc T^cc 3I b b 



3X^ C — 3X; 



6c 



(B18) 



The integration of in the soft region with dimensional regularization can be found in 
Ref. 54]. and with massive gluon can be found in Ref.[55j] 
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APPENDIX C: THE TOTAL DECAY WIDTH 



The decay width at NLO in a s is 

FnLO = r_LQ + r Virtual + ^ Real- 



(CI) 



The LO decay width of the color singlet piece has been given in Eq.([2]). The D-dimension 
LO decay width is 

r^TOX.) - = ^mzllfl f ( IV Y <T| (3 5l , 1 )|T> , (C2) 

81m^ 2r[| — ej \P m b / 



where r = m c /m b and /3 = Vl — ?" 2 - 

The rV ir4tta / of the color singlet piece can be written as 



^Virtual{ 3 S^l) Q> s V 2 + 2 



32i7r 2 



it r z + V — 2 I \e 



\m b m c ) ) 



+ 



16m 2 



9 (2m 2 + m c 2 ) 



3 (2m 2 + ml) 9 

A (m 6 ) ( 4 - ^| ) - 3A (m c ) - 2 (3£? (p c , 0, m c ) (4m 2 + m 2 ) + 
\ m^ / \ 



52 

+ 



.Bo (2p b , m b , m b ) (2m 2 , + m c) _ 35 o (2p&, m c , m c ) (3m 2 + m 2 ) + 3B (Pb, 0, m b ) 
+12 ({D - 2)m 2 + m 2 J (c (p b , -p b , 0, m b , m b ) m 2 + 

Co (Pc ~Pc, 0, m c , m c ) (2m 2 - m 2 ) 

For the real corrections, we should introduce the Mandelstam variables 



///; 



(C3) 



In the rest frame of T, 



S23 = (Pc+Pc) 

s 3 4 = (Pc + kf 



s 2 3 = 4ml ~ 4m b fc c 



(C4) 



(C5) 



s 34 = Ami + m l- 4m b p®. 
The contribution of real corrections for the color singlet piece is 

diW 3 Si,i) lejTjg^QjTj^ x 

dfc°dp° 243m^ (m 2 — s 34 ) 2 (—4m 2 — m 2 + s 23 + s 34 ) 2 

2m® - 8s 34 m^ + (3s 2 3 + 4s 34 ,s 23 + 12s 2 4 ) - (s^ + 2s 34 s 2 3 + 8s 34 s 23 + 8s 34 ) m\ + 
64m^ (3m 2 - s 34 ) + s 34 (s 23 + s 34 ) (s 2 3 + 2s 34 s 23 + 2s 34 ) + 



16m t 



lm\ - (5s 23 + 6s 34 ) m 2 + s 34 (s 23 + 3s 34 ) 



+ 



Ami 



Am\ - 4 (2s 23 + 3s 34 ) m\ + (3s 2 23 + 4s 34 s 23 + 12s^ 4 ) m 2 - s 34 (s 23 + 2s 34 ) 2 \ . (C6) 
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The NLO decay width of color singlet piece is 



J- NLO 



ana 
81m 2 



;T|0( s ^)|T>|vT^^l + y) 



(32 - 8r 4 )Li 2 (x /3 ) + (16 - 4r 4 ) {lA 2 {-xp) + ln(^) ln(l - x p ) 

/ \ 9r 2 

+ (2 + r 2 )Vl-r 2 (6 \n{x p ) - 8 ln(l - x p ) - 4 ln(l + X/j) J + (3 + — ) VT 

7r 4 

+ (-12 + 2r 2 + Info) + (8 - 2r 4 ) Info) ln(l + a^) 



(C7) 



where /3 = Vl 3 ^ 2 and x p = (l- f3)/{l +J) = (1 - >/l^r 5 )/(l + Vl 3 ^). The ~^C F 
term is due to the QCD correction to Y[&&( 3 S' ljl )] — > 7*, while the j*-Cf[- ■ ■} term is due to 
the QCD correction to 7* — » cc. This is the same as the known next-to-leading order result 
of e+e~ -> 7* -> cc (HEe3- 

The LO decay width of the color octet piece has been given in Eq. (|2ip . The D-dimension 
LO decay width is 



r LO [T( 3 5 1 , 8 ) 



cc 



a 2 s (3 (D-2 + r 2 



71 



6m 2 2T[ 
The Tyirtuai of the color singlet piece can be written as 
r V irtuai( 3 S 1>8 ) r 2 + 2 I 117/1 



47T/X 

/3 2 m 2 



2 \ e 



<T|G( 3 5 li8 )|T). (C8) 



r LO ( 3 ^i 



9vr r 2 + D-2 



- 7b + In (Air) - 18 log 



mlm 2 c 



+ 



202m 2 - 193m 2 



+ 



2m 2 



A (m b ) (35m 2 - 2m 2 ) - 12A Q (m c ) (6m 2 + m 2 ) 



4 (2m 2 + m 2 ) ' m 2 (2m 2 , + m 2 ) - 
15B (pt, 0, m 6 ) (7m 4 + 12m 2 m 2 ) - 905 {2p b , 0, 0) m 2 (2m 2 + m 2 ) + 
4S (2p b , m 6 , m b ) m 2 b (2m 2 + m 2 ) + 12B (2p b , m c , m c ) (m 4 + 3m 2 m 2 + m 4 ) - 
96^0 (Pc 0, m c ) m 2 (4m 2 + m 2 ) - IO8E0 (Pb ~ Pc, m b , m c ) m 4 - 
324C {pb, ~Pb, m b , 0, 0) m\ (2m 2 + m 2 ) - 216C {Pc ~Pc, m c , 0, 0) m\ (2m 2 + m 2 ) 
24C (Pb, -Pi, 0, m 6 , m 6 ) m 4 ((£) - 2)m 2 + m 2 ) - 
216C (Pc, Pfo, 0, m c , m 6 ) m 4 ((£> - 2)m 2 b + m 2 ) - 
24C (p c , ~Pc, 0, m c , m c ) m 2 , (2m 2 ; - m 2 ) ((£> - 2)m 2 b + m 2 ) - 

432D (p b , p b - p c , -p b , m b , 0, m c , 0)ml((D - 2)m 2 + m 2 ) 



(C9) 
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The contribution of real corrections for the color octet piece is 

(T|0( 3 S 1)8 )|TK 



dk°dp° c 



18m 4 (s 23 - 4m^) 2 



S34) 2 (-4m 2 - m 2 + s 23 + S34) 2 



x 



64m 4 + 4 (9m 2 - 8s 23 - 9s 34 ) m 2 + 9m 4 + As\ 3 + 9s 2 3i + 9s 23 s 34 - 9m 2 (s 23 + 2s 34 ) 



x 



2m® - Bs 34 m 6 c + (3s 2 23 + 4s 34 s 23 + 12s 34 ) m 4 - (s\ 3 + 2s 34 s 2 23 + 8s 34 s 23 + 8s 3 u ) m 2 + 
64m^ (3m 2 - s 34 ) + s 34 (s 23 + s 34 ) (s 2 3 + 2s 34 s 23 + 2s 34 ) + 



lQm h 



7m 4 - (5s 23 + 6s 34 ) m 2 + s 34 (s 23 + 3s 34 ) 



+ 



Am 2 b 



Ami ~ 4 (2s 23 + 3s 34 ) K + (3s 2 3 + 4s 34 s 23 + 12s 2 4 ) m 2 - s 34 (s 23 + 2s 34 ) 2 \ . (CIO) 



The NLO decay width of color octet piece is 



ag(T|Q( 3 g 1|8 )|T) 
216m 2 



bml 

2tt 2 /3 4 - 478/? 3 + lAir 2 /3 2 + 1614/3 - 60tt 2 - 150/3 (f3 2 - 3) In + 
6 (/3 2 - 3) ln(2) [90/3 + (/3 2 - 17) ln(2)] + 408/3 (/3 2 - 3) ln(/3) - 



12/3 5 + 



-6/3 (35/3 2 - 87) log(l - (3 2 ) - 6 (/3 6 - 13/3 4 + 4/3 2 + 60) log 



1-/3 
/3 + 1 



+ 



61og(/3 + l) log(2)/3 4 -51og(16)/3 2 + (-2/3 4 + 31/3 2 - 75) log(/3 + 1) + 51 log(2) + 
6 log(l - /3) [ - log(8)/3 4 + 60 log(2)/3 2 + 3 (/3 4 - 20/3 2 + 51) log(/3 + 1) - 153 log(2) 



12 (/3 2 - 3) 9Li 2 (/3) + 9Li 2 



/3 + i; + ^ 2 - 17 ) ( 2Ll2 V/3 + i 



2/3 



Li 2 



/3 + 1 



12 (/3 4 - 20/3 2 + 51) log(/3) log 



1-/3 
/3+1 



(Cll) 
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